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The relativistic local density approximation (LDA) for the Coulomb exchange functional in nuclear 
systems is presented. This approximation is composed of the well-known Slater approximation in the 
non-relativistic scheme and the corrections due to the relativistic effects. Its validity in finite nuclei 
is examined by comparing with the exact treatment of the Coulomb exchange term in the relativistic 
Hartree-Fock-Bogoliubov theory. The relativistic effects are found to be important and the exact 
Coulomb exchange energies can be reproduced by the relativistic LDA within 5% demonstrated by 
the semi-magic Ca, Ni, Zr, Sn, and Pb isotopes from proton drip line to neutron drip line. 
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I. INTRODUCTION 

The Coulomb interaction between protons is one of the 
most important building blocks in atomic nuclei P, . It 
is the main source of the charge symmetry-breaking in 
the nuclear Hamiltonian. Globally, the Coulomb inter- 
action is the driving force for determining the /3-stability 
valley away from the N = Z line. Specifically, the 
Coulomb interaction plays a crucial role in understand- 
ing the Coulomb displacement energies, isospin mixing, 
proton emission, fission barriers, a-decay energy release 
Q a in superheavy elements, and more. Therefore, the 
nuclear density functional theory (DFT) that models the 
effective strong interactions has to be accompanied by an 
explicit functional form for the Coulomb interaction 0). 

Generally speaking, the Coulomb interaction is the 
best known part of the nuclear Hamiltonian, and both 
its direct (Hartrce) and exchange (Fock) terms can be 
exactly calculated at the mean-field level. However, the 
exchange term is very cumbersome to include due to 
the non-locality of the corresponding mean field, espe- 
cially for the deformed nuclei. In the spirit of the DFT 
of Kohn and Sham [|| , the local density approximation 
(LDA) has been intensively adopted in the calculations 
of finite nuclei since the 1970s Q. So far, in almost all 
of the non-relativistic (NR) self-consistent Hartree-Fock 
calculations, the Coulomb exchange energy and the corre- 
sponding single-particle potential are evaluated within a 
local scheme by using the so-called Slater approximation 
Q . The validityof this approximation has been discussed 
in the Skyrme |3d^| and Gogny [Io| approaches. 

During the past decades, the nuclear covariant den- 
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sity functional theory (CDFT) at the Hartree level, also 
known as the relativistic Hartree (RH) or relativistic 
mean-field (RMF) theory, has received much attention 
due to its successful description of many nuclear phe- 
nomena The Lorentz covariant form of the the- 
ory itself guarantees the self-consistent spin-orbit poten- 
tial and also puts stringent restrictions on the number 
of parameters in the corresponding functionals without 
reducing the quality of the agreement with experimental 
data. In order to retain the required simplicity in the 
RH framework, the non-local Coulomb exchange term 
is usually neglected. Its effects on the binding energies, 
radii, etc., are assumed to be absorbed into the effective 
coupling strengths through the phcnomcnological fit of 
the model. However, this prescription of neglecting the 
Coulomb exchange term is not always valid. One exam- 
ple is the isospin symmetry-breaking corrections to the 
superallowed j5 decays [l8|], which are crucial for testing 
the unitarity of the Cabibbo-Kobayashi-Maskawa matrix 

Recently, the density-dependent relativistic Hartree- 
Fock (RHF) theory [20, [2l| has been developed and 
achieved equivalent success in describing the ground- and 
excited-state properties as the RH theory. It has been 
shown that the meson exchange terms play very impor- 
tant roles in the nucleon effective mass splitting [20J, 
symmetry energies [22| , spin and pseudospin symmetries 
[23l |24j . shell structure and its evolutions [2ll I25I - 
deformation [2{|, and spin- isospin resonances [30l 
Subsequently, a unified and self-consistent description 
of both RHF mean field and pairing correlations, i.e., 
the relativistic Hartree-Fock-Bogoliubov (RHFB) theory 
[HI, [33| , has been achieved for the exotic nuclei far from 
the /3-stability valley. In the frameworks of RHF and 
RHFB, for the first time, the non-local Coulomb ex- 
change term has been taken into account exactly, which 
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is crucial for understan ding the isospin corrections to the 
superallowed /3 decays ^8j . 

With the simplicity of the local CDFT and the success 
in treating the Coulomb terms exactly, it is worthwhile 
to explore the relativistic LDA for the Coulomb exchange 
term, a topic less known in the nuclear physics commu- 
nity. The relativistic LDA may be a simple approach 
to implement the missing Coulomb exchange effects in 
the RH theory, yet keeping the merits of locality. The 
attempt in this direct is further supported by a recent 
success in folding the meson exchange terms into a local 
equivalent scheme [HI]. 

The relativistic LDA has been explored for electronic 
systems for the past decades. The extension of the Ho- 
henberg and Kohn theorem 13511 to relativistic systems 
was first formulated in Ref. j36j by utilizing a quan- 
tum electrodynamics (QED)-bascd Hamiltonian with the 
four-current. Within the no-sea approximation, which 
is also called no-pair approximation in atomic physics 
for neglecting all effects due to the creation of particle- 
antiparticle pairs, the total energy of the system can be 
expressed as a functional with respect to the four-current 
ji*(r) = (p(r)J(r)) [II [H, instead of only p(r) in the 
NR case. In complete analogy with the NR case, the rela- 
tivistic LDA for the Coulomb exchange energy functional 
Ecexlj^i?)] is constructed by using the model of relativis- 
tic homogeneous electron gas. As the space-like compo- 
nent j vanishes in the homogeneous systems, the full j^- 
dependence of the exact exchange energy functional is 
reduced to a pure density-dependence in the relativistic 
LDA. The relativistic corrections, i.e., the differences be- 
tween E^ A [p(r)} and E^ BA [p(r)}, are shown to be 
substantial, see Ref. [3t| for details. 

In this paper, the relativistic LDA for the Coulomb in- 
teraction in nuclear physics will be presented. The main 
focus will be its validity by comparing with the exact 
RHFB results and the relativistic corrections to the tra- 
ditional Slater approximation. In Sec. [TTl the Coulomb 
energies in DFT and the corresponding LDA for the ex- 
change terms will be presented. The numerical details 
and the results will be presented and discussed in Sec. lIIII 
Finally, a summary will be given in Sec. IIVI 



is simply a functional of p p (r), i.e., 
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and the corresponding Hartrec potential for protons 
reads 



Vbdir(r) = — 7-t = e 
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In contrast, the exchange term of Coulomb energy, 

,^(r)^(r)^(r')V,(r') 
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r - r 



(4) 

is non-local and thus more complicated to calculate. 

In the spirit of the LDA, the exchange energy density 
of an inhomogencous system with density p(r) is locally 
approximated by the exchange energy density of a homo- 
geneous system with density n = p(r). In the homoge- 
neous nuclear matter, the single-particle wave functions 
are the plane- wave solutions of the Schrodinger equation, 
so that the Coulomb exchange energy shown in Eq. (|4]) 
per unit volume can be calculated analytically as 
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where n p is the proton density. Therefore, by using the 
LDA, the Coulomb exchange energy in finite nuclei is 
determined as 
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This is the well-known Slater approximation [6[ widely 
used in the non-relativistic DFT. The corresponding 
single-particle potential for protons reads 



t/LDA / \ _ SE^ A [ Pp (r)] 

V Ce X [ r — 1 73 



1/3 



?W 3 (r). (7) 



II. COULOMB ENERGIES AND LOCAL 
DENSITY APPROXIMATIONS 

Before introducing the relativistic LDA, it is illuminat- 
ing to briefly recall the Coulomb energies and the corre- 
sponding LDA for the exchange term in the NR systems. 

In the non-relativistic scheme, the proton density dis- 
tribution in nuclei reads 



PpO) 



(1) 



with the single-particle wave functions ipi (r) and occupa- 
tion probabilities vf . The direct term of Coulomb energy 



In the relativistic scheme, the variables of the energy 
functional include not only the density distributions but 
also the currents. The proton density distribution and 
currents read 

Pp (v) = ^^,(rhVW, (8a) 

i 

j p (r) = ]T^(r)7<Mr), (8b) 

i 

respectively. It should be emphasized that the summa- 
tions over i include only the particles in the Fermi sea, 
i.e., the no-sea approximation, which is the so-called no- 
pair approximation in atomic physics [39j . 
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Correspondingly, the direct term of Coulomb energy 
reads 



^Cdir 
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where the first and second terms are the contributions 
from the time-like and space-like components of the 
Coulomb field, respectively. They are called Coulomb 
and transverse contributions in Ref . [39| . For the systems 
with time-reversal symmetry, the contribution from the 
space-like component vanishes. Then, Eq. © has the 
same structure as Eq. ([2]). 

On the other hand, the non-local exchange term of 
Coulomb energy is j39| 



e 2 <A 



d 3 rd 3 r 



,cos(|£i - £j||r - r'| 



x^(r)7%-(r)^(r')7^(r'), 



(10) 



where £i are the single-particle energies. 

In the relativistic homogeneous nuclear matter, the 
single-particle wave functions are the plane-wave solu- 
tions of the Dirac equation. As a result, the time-like 
component e£ 0X and space-like component e§ ex of the 
Coulomb exchange energy shown in Eq. (| 10p per unit 
volume can be respectively expressed as [2a, |3_9j 



egex = e C ex$(/?) and eg e 



e C cx$(/3), 
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by the NR one ec ex in Eq. ([5]) together with 
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where /? = (Stt 2 ^) 1 / 3 /M ', i] = </l + /3 2 , and M is the 
proton mass. The details can be found in Ref. [1^ and 
the references therein. 

From Eqs. ([TT|) and ([T2"j) . the relativistic results are 
identical to the NR counterparts at the zero density limit. 
The relativistic corrections to the Coulomb exchange en- 
ergy increase with the density. For the symmetric nuclear 
matter, n p w 0.08 fm~ 3 at the saturation density, one has 



P ~ 0.28, (3 2 ~ 0.078, and /? 4 - 0.006. Therefore, the 
relativistic corrections to the Coulomb exchange energy 
are expected to be substantial in nuclear systems, and 
the contributions from /3 4 and higher-order terms can be 
neglected. Accordingly, the relativistic effect from the 
space-like component is 5 times of that from the time- 
like component, as shown in Eq. (|12[) . 

From Eq. (fTTj) . up to the order of (3 2 , the Coulomb 
exchange energy in the relativistic LDA is 
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The corresponding contribution to the single- particle po- 
tential for protons reads 



(|) 1/3e V 3 w + (^) 



(14) 



III. RESULTS AND DISCUSSION 



In the following, the properties of finite nuclei are cal- 
culated in the density-dependent RHFB theory [32[ with 
PKA1 [2l| and D1S [1| effective interactions for the 
particle-hole and particle-particle channels, respectively. 
The RHFB equations are solved on a Dirac Woods-Saxon 
basis [4l| within a spherical box of radius i? max = 20 fm, 
and the numbers of positive and negative energy levels 
for each (Z,j)-state are fixed to Np = 28 and Np, = 12, 
respectively. 

In the standard RHFB calculations, the Coulomb ex- 
change energy and the corresponding non-local single- 
particle potential are calculated exactly. They are la- 
beled with the superscript "exact". The correspond- 
ing results obtained by using the relativistic LDA shown 
in Eqs. (TIB")) and (| L4|) are labeled with the superscript 
RLDA. In order to quantitatively investigate the rela- 
tivistic corrections in Eq. (fTTj). the results of the non- 
rclativistic LDA shown in Eqs. (JB]) and ([7]) labeled with 
the superscript NRLDA are also presented for compari- 
son. The relative deviations of the approximate Coulomb 
exchange energies from the RHFB results are defined as 



rpLDA pexact 

AEr = Ccx Ccx 
Ccx pcxact 
^Ccx 



(15) 



In order to exclude the effects due to the self- 
consistency, one-step calculations have been performed 
to investigate the effects of the LDA approximation on 
Coulomb exchange energies, i.e., E^l^ A and E^^ BA arc 
respectively obtained with the relativistic and NR LDA 
by the same converged proton density distributions p p (r) 
given by the self-consistent RHFB calculations. By tak- 
ing the even-even Pb isotopes from proton drip line to 
neutron drip line as examples, the Coulomb exchange 
energies Eq qx arc shown as a function of mass number 
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FIG. 1: (Color online) Coulomb exchange energies in Pb iso- 
topes calculated by RHFB theory [H with PKA1 [2l|. The 
results obtained with the same proton density distributions 
p p (r) but within the relativistic and non-relativistic (NR) lo- 
cal density approximations (LDA) and their relative devia- 
tions from the exact results are shown in the upper and lower 
panels, respectively. 
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FIG. 2: (Color online) Proton density distributions p p (r) in 
188 Pb, 208 Pb, 228 Pb, and 248 Pb by RHFB theory with PKA1. 



A in the upper panel of Fig. [T] The corresponding rela- 
tive deviations AEq cx dchned in Eq. (|15|) are shown in 
the lower panel. In addition, the proton density distribu- 
tions p p (r) in the nuclei 188 Pb, 208 Pb, 228 Pb, and 248 Pb 
are illustrated in Fig. [2J 

In general, the magnitudes of the Coulomb exchange 
energies i?ccx decrease with increasing mass number, i.e., 
with increasing size. For each nucleus, the magnitude of 
Eqcx is overestimated by the NR LDA or the traditional 
Slater approximation, and substantially improved when 
the relativistic corrections are taken into account. From 
the relative deviations (see the lower panel of Fig. [l| , the 
relativistic correction improvements range from 4.5% in 
182 Pb to 3.5% in 266 Pb. The isotopic dependence of the 
relativistic effect is due to the proton density distribu- 
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FIG. 3: (Color online) Proton single-particle energy shifts 
caused by the RLDA in 208 Pb. The Fermi energy is shown as 
the vertical dash-dotted line. 



tions as shown in Eq. (|13[) . in particular for the interior 
region of nuclei. As shown in Fig. [51 the inner proton 
density decreases with the neutron number, which re- 
sults in a smaller relativistic correction. Accordingly, the 
relative deviations AEcex due to the relativistic LDA de- 
crease from 3.6% in 182 Pb to 1.7% in 266 Pb. As shown 
in Fig. [2] the inner proton density in the heavy nuclei 
is almost flat and forms a near uniform density region. 
This is the reason why the LDA works well. 

In order to investigate the effects due to the self- 
consistency, the self-consistent RLDA calculations have 
been performed. In Fig. [3l the proton single-particle en- 
ergy shifts e^ LDA — e° xact in 208 Pb are shown as a func- 
tion of the single-particle energy £p Xact . It is seen that, in 
the relativistic LDA, the occupied and unoccupied pro- 
ton states are pushed up and down, respectively. The 
energy shifts are around +400 keV for the deeply bound 
states, +100 keV for the states near the Fermi surface, 
and —200 keV for the unoccupied bound states. A very 
similar shifting behavior was shown in Ref. 0] in the 
NR Skyrme Hartree-Fock framework. As the occupied 
states are pushed up due to the RLDA, the total kinetic 
energy of protons increases by 2.73 MeV, which is com- 
pletely cancelled out by the self-consistent effects of the 
neutron kinetic (—2.12 MeV), mesons (—1.18 MeV), and 
Coulomb direct (+0.58 MeV) energies. The same kind 
of self-consistent cancellation between the kinetic energy 
and other channels was also noted in Ref. @. As a result, 
the total energy difference between the RLDA and RHFB 
calculations, —0.97 MeV, is mainly due to the difference 
in the Coulomb exchange energies, —0.99 MeV. 

In Fig. [H we show the systematics of calculated re- 
sults using the self-consistent relativistic LDA for the 
semi-magic Ca, Ni, Zr, Sn, and Pb isotopes from proton 
drip line to neutron drip line. The traditional doubly 
magic nuclei are marked by the open symbols. Com- 
paring with the exact Coulomb exchange energies calcu- 
lated in RHFB, it is found that the relative deviations 
A-Ecex introduced by the relativistic LDA are less than 
5% for all these five semi-magic isotopes. In addition, 



5 




-5 I I i I i I i I i I i I i 1 

40 80 120 160 200 240 280 

A 

FIG. 4: (Color online) Relative deviations of the Coulomb 
exchange energies by the self-consistent RLDA (solid lines) 
and NRLDA (dashed lines) calculations for Ca, Ni, Zr, Sn, 
and Pb isotopes. The traditional doubly magic nuclei are 
denoted as open symbols. 

from the differences between the solid and dashed lines, 
the relativistic corrections to the LDA are found to play 
substantial roles in improving the results by 3 ~ 5%. 

Therefore, one can conclude that the relativistic ver- 
sion of LDA shown in Eqs. (fl~3"|) and (JUJ) for the Coulomb 
exchange term in nuclear CDFT is very robust and 
promising. In particular, the relativistic corrections to 
the traditional Slater approximation shown in Eqs. 
and (|12p are very important. 

IV. SUMMARY 

In this work, the relativistic LDA for the Coulomb ex- 
change functional in nuclear systems is presented. This 



approximation is composed of the well-known Slater ap- 
proximation in the NR scheme and the corrections due 
to the relativistic effects. The validity of the relativistic 
LDA in finite nuclei calculations is examined by com- 
paring with the results of the RHFB theory, where the 
non-local Coulomb exchange term is treated exactly. It 
is found that the relative deviations of the Coulomb ex- 
change energies in the self-consistent RLDA calculations 
are in general less than 5% for semi-magic Ca, Ni, Zr, 
Sn, and Pb isotopes from proton drip line to neutron 
drip line. For the proton single-particle energy shifts, 
the relativistic LDA pushes the occupied and unoccu- 
pied states upward and downward, respectively, which is 
in agreement with the previous NR Hartree-Fock calcula- 
tions. Finally, it is also worthwhile to emphasize that the 
relativistic corrections to the LDA are found to play sub- 
stantial roles in improving the agreement with the exact 
results by 3 ~ 5%. 
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